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Abstract. We find a transformation which relates a new third-order inte- 
grable nonlinear evolution equation, introduced recently by Qiao, with the 
well-known modified Korteweg-de Vries equation. Then we use this transfor- 
mation to derive smooth soliton solutions of the new equation from the known 
rational and soliton solutions of the old one. 



1. Introduction 

In the present paper, we study the following third-order integrable nonlinear 
evolution equation: 

This equation was introduced recently by Qiao [1] as the second positive member 
in a new completely integrable hierarchy. The Qiao equation — we will refer to ([1]) 
in this way — possesses a Lax representation and a bi-Hamiltonian structure [TJ [2] ■ 
Moreover, the Qiao equation may be reduced from the two-dimensional Euler equa- 
tion by an approximation procedure, and its solutions may be useful to construct 
new solvable potentials in Newtonian dynamics and to model electrophysiological 
phenomena in neuroscience pQ. 

There is an interesting problem concerning the soliton solution of the Qiao equa- 
tion. In [TJ where the traveling wave solutions of (JJ were studied, it was stated 
that this equation has no smooth solitons, and only piecewise smooth, but not 
smooth, continuous traveling wave solutions were found in those papers. On the 
other hand, in [3], where the bifurcations of traveling wave solutions of (Q]) were 
studied, it was stated that this equation does possess smooth soliton solutions, and 
some of those solitons were shown in Figure 5 of that paper. We believe that the 
complexity of investigation of the Qiao equation's solitons is caused mainly by the 
fact that they can be expressed only in a parametric form. In the present paper, we 
also address the problem on smooth soliton solutions of ((TJ, but we do not use the 
traveling wave reduction to an ordinary differential equation. Instead, we follow 
the way used in [3] for the short pulse equation: we transform the studied new 
equation into a known old one, and then transform known explicit expressions for 
solutions of the old equation into sought parametric expressions for solutions of the 
new one. 

The paper is organized as follows. In Section [2] we find a transformation which 
relates the Qiao equation with the well-known modified Korteweg-de Vries (mKdV) 
equation. Then, in Sections |3] and |4] we apply this transformation, respectively, to 
the rational solution and to the soliton solution of the mKdV equation, and in this 

1 



2 



S. SAKOVICH 



way we obtain three types of smooth solitons of the Qiao equation as well as three 
types of its smooth antisolitons. Section [5] contains some concluding remarks. 

2. Transformation to the mKdV equation 

Let us try to transform the new Qiao equation ([!} into one of well-studied old 
equations. We do this in the way used in [5] for the Fujimoto-Watanabe equations; 
see also [5J [7] for more theory and applications of the transformations involved. 

First of all, we transform the nonlinear evolution equation ([T]), which has the 
form m t = —rn~ 3 m xxx + (lower-order terms), into an equation of the form Ut — 
u 3 u xxx + (lower-order terms). To do this, we change in ([l} the dependent variable 
m as 

m{x > t) = ^jy (2) 

the result being the polynomial evolution equation 

u t = v?u xxx + 3u 2 u x u xx — u 3 u x . (3) 

This equation is just a special case of the item (II. 2) in the Fujimoto-Watanabe 
classification [8], Ut = u 3 u xxx + 3u 2 u x u xx + 4au 3 u x with the arbitrary constant 
a fixed to be —1/4. Therefore, according to [5], it is possible to relate the Qiao 
equation, through its previously known form ([3]) , with the Korteweg-de Vries equa- 
tion. We will see soon, however, that it is more convenient to transform the Qiao 
equation into the mKdV equation, in order to easily derive solutions of the former 
from solutions of the latter. 

Next, we take into account that the evolution equation (j3j can be written as 

u t = u 2 (uu xx + u 2 x - iu 2 ^ , (4) 

which means, according to (51 [5], that the transformation 

x = v(y, t), u(x, t) = v y (y, t) (5) 
relates with the third-order evolution equation 

V t = Vyyy ~ ^ . (6) 

This point should be explained in more detail. It is easy to see that any function 
v(y,t) satisfying ([6]) determines via ([5j a function u(x,t) satisfying In the 
opposite direction, when we substitute the expressions ([5]) to ([3]), we obtain a 
fourth-order equation for v(y, t), which, owing to the form of Q, can be written as 

Hl^ + l v A =0 (7) 
and then can be integrated over y, the result being the third-order equation 

V t = Vyyy ~ + f (t)Vy , (8) 

where f(t) is an arbitrary function. This means that any solution u(x,i) of ([3]) is 
represented parametrically, via ([SJ with y being the parameter, by a solution v(y, t) 
of (|SJ) with some function fit). However, we can always fix this function f(t) to be 
zero, by the corresponding change of the parameter y, y i— > y + J f(t)dt, which has 
no effect on the represented solution u(x,t). Consequently, all solutions of ((3j are 
represented by all solutions of ((6|) via the relations ([5]). 
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Finally, we introduce the new dependent variable w, 

w(y,t) = v y (y,t), (9) 

which is convenient because the relation 

u{x,t)=w(y,t) (10) 

holds due to §5$. According to © and ©, the function w(y,t) is any solution of 
the mKdV equation 

3 

W t =Wyyy ~ ~W Wy, (ll) 

whereas the function v(y,t) is determined by w(y,t) via the relations v y = w and 
v t = Wyy — ^w 3 up to an arbitrary additive constant of integration. Of course, there 
is no need to further transform (|11|) into the Korteweg-de Vries equation, because 
it would only complicate the already found way to solve the Qiao equation. 

We can summarize the obtained transformation, which relates the equations 
(TTJ) and (ITT]) to each other, as follows. All solutions of the Qiao equation (TJJ arc 
represented parametrically by all solutions of the mKdV equation (jlip via the 
relations 

mix, t)= 7 . , (12) 

x = v(y,t): { i 3 (I.T) 

where the independent variable y serves as the parameter. Of course, this transfor- 
mation is local, in the sense that it maps a local solution of (fTTj) to a local solution 
of (TTJ), and for every particular w(y, t) some extra analysis is required to see whether 
the corresponding m(x,t) is real- valued, smooth, and global. 

In the next two sections, in order to show how the transformation flTZt (fl~3|) 
works, we apply it to two simplest nontrivial solutions of the mKdV equation (|llj) : 
the rational solution 

4a 3 2 

w{y, t) = a + -5-5 — -, r = y--at, (14 
a z r z — 1 2 

where a is an arbitrary nonzero constant, and the soliton solution 

2k 2 /3 




w(y,t) = a + =^——— , s = y- I ^ra 2 - k 2 ) t, (15) 

±Vfl — k ^ cosh(fcs) — a \4 J 

where a and k are arbitrary constants, k ^ and a 2 — k 2 ^ 0. One can notice that 
the expressions (|14[) and (|15p are somewhat different from the ones well known in 
the literature [5]. The reason is that the mKdV equation is usually studied in its 
focusing form, with a positive coefficient at the nonlinear term, whereas we have the 
defocusing case in (fTTj) . and the expressions (|14j) and (|15|) can be obtained from the 
well-known ones through a complex-valued scale transformation of the dependent 
variable and a change of notations for the constants. 
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3. Transformation of the rational solution 

Let us derive a solution of the Qiao equation ([l} from the rational solution (|14[) of 
the mKdV equation (TTT1) by means of the transformation (fl ^ - (fT"3"|) . We substitute 
the function w(y, t) given by (fl4]) to the relations (fT2j) and (|T3|) . make the integration 
required in (|13|) . and fix the arbitrary additive constant of integration so that to 
have v = for j/ = t = 0. In this way, we obtain the following expressions: 

1 — a 2 r 2 

m ( x ^)=o — —3-^"' ( 16 ) 
x = a 3 t + ar - 4tanh -1 (ar), (17) 

where tanh -1 stands for the hyperbolic arctangent. These expressions contain y via 
the variable r — y — ^a 2 t only, therefore we can use r instead of y as the parameter, 
for simplicity. The function m(x, t), determined parametrically by (116l) - (|17[) with 
r being the parameter, does satisfy (JXJ) , as one can prove by direct computations. 

Until now, we worked with the solution (p~4]t in a purely algebraic way, and 
the quantities a, y, t and w could be even complex-valued there. However, in order 
that the obtained solution m(x, t) (fTB| - (fT7| be a real-valued function of real-valued 
variables, the quantities a, y and t must be real- valued, and the condition 

- 1< ar < 1 (18) 

must be satisfied, which is caused by the presence of tanh -1 in (TT7)) . The domain 
of variation of x is — 00 < x < 00 due to (flTf and (fl8|) . The function m(x,t), due 
to (|16| and (fi"8|) . is everywhere smooth; it is positive, with one maximum equal to 
l/(3a), when a > 0; it is negative, with one minimum equal to l/(3a), when a < 0; 
and its asymptotes at x — > ±00 are zero. 

We can see easily that the obtained solution ([TO]) — (TT5T) is a traveling wave solu- 
tion, m(x, t) = F (a; — a 3 <) with the function F being determined parametrically. 
The shape of this solitary wave is close to the shape of a hyperbolic secant curve 
but not identical to it, as shown in Figure [TJ The speed of this solitary wave is 
equal to a 3 ; solitons (positive waves with a > 0) move to the right, and antisoli- 
tons (negative waves with a < 0) move to the left, as shown in Figures [5] and O 
respectively. This antisymmetry m t— > — m, x <— ¥ —x between positive and negative 
solutions of the Qiao equation ([1]) corresponds, via the transformation (fT^ — (fl3|) . 
to the symmetry w >->• — w of the mKdV equation (|11| . 

To conclude the section, we point out the following two discovered phenomena. 
First, a rational solution of one nonlinear equation is transformed into a smooth 
soliton solution of another nonlinear equation. Second, only a part of the rational 
solution (fl4)) of the mKdV equation — the part bounded by the condition (fT8|) — is 
used to produce the entire smooth real solution of the Qiao equation. 



4. Transformation of the soliton solution 

Now, let us proceed to transformation of the soliton solution (TT5)) of the mKdV 
equation (flT]) . The cases of "+" and "— " signs in (fl5)) will be considered separately. 

We substitute the function w(y, t) (fT5|) . with the choice of "+" sign in it, to the 
relations (fT2|) and (fT"3"|) . make the integration required in (fT3| . fix the constant of 
integration so that to have v = for y = t = 0, and obtain in this way the following 




FIGURE 3. The moving antisoliton P^ ]) -([I5 ]) with a = -1: t = 
(solid) and t = 10 (dashed). 
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expressions: 



a — V a 2 — fc 2 cosh(fcs) 

m(x,t) = - v ; , (19) 

2k 2 — a 2 + a\J a 2 — k 2 cosh(fcs) 

a 2 -k 2 )t + as - 4 tanh -1 f -j- (a + \/a 2 - fc 2 ) tanh |) . (20) 



x = a i 

These expressions contain y via the variable s — y — (§a 2 — fc 2 ) i only. Therefore we 
can use s instead of y as the parameter, for simplicity, when we consider (|T9l) - ([20|) 
as a parametrically determined solution m(x,t) of the Qiao equation ([TJ. This is a 
traveling wave solution, m(x, i) = G (x — a (a 2 — fc 2 ) i) with the function G being 
determined parametrically, and the speed of the wave is equal to a (a 2 — fc 2 ) . Of 
course, in order that the obtained solution be a real solution, the quantities a, fc, s 
and t must be real- valued in the right-hand sides of (|T9")) and (j2U)) (we take fc > 
without loss of generality), and the condition fc < \a\ must be satisfied due to the 
presence of the square root. Further constraints on s and fc, caused by the presence 
of tanh -1 in (|20|) and the denominator in (fl"9|) . depend on the sign of a. 
When a > 0, the parameter s must lie in the interval 

|s| < j tanh" 1 k < fc < a, (21) 

fc a + Va 2 - fc 2 

in order that x (|20l) be real- valued. With this constraint on s, the domain of varia- 
tion of x is — co < x < co, and the function m(x, t) (|19|) is everywhere smooth and 
positive, with one maximum equal to 1/ (a + 2-\/ a 2 — fc 2 ) and zero asymptotes at 
x — > ±co. The soliton, determined by (fT9"]) - (l2"T1) . moves with the speed a (a 2 — fc 2 ) 
to the right, as shown in Figure SJ When fc — > 0, the height and the speed of this 
soliton tend to l/(3a) and a 3 , respectively. In fact, the soliton (TTB|) ~([T5 |) with a > 0, 
derived from the rational solution of the mKdV equation in Section [3l is a limiting 
case of the soliton p9|) -([2~T |) with fc — > 0, as one can prove easily. On the other hand, 
when fc — > a, the height and the speed of the soliton (fT9"]) - (|2"Tj) tend to 1/a and 0, 
respectively. Figure [S] shows how the shape of the obtained soliton depends on the 
ratio k/a. We see that the soliton's top widens and flattens, when fc is very close 
to a, that is when the soliton's motion is extremely slow. However, the soliton's 
shape changes very little, when fc runs over the wide interval < fc < 0.97a. This is 
shown in Figure [6J where we have adjusted the values of a so that to have solitary 
waves of equal heights. It is interesting that, of the two solitary waves shown there, 
the one which is only a little bit wider moves more than two times slower than the 
other one. 

When a < in (| 19|) and (|20)) . the variable x ([20]) is real- valued on the whole 
interval — co < s < co for any fc, < fc < —a, and the domain of variation of 
x is —co < x < co, but the function m(x,t) (I19p has one or two singularities if 
fc = — \/3a/2 or fc > — \/3a/2, respectively. The expressions (fT9]) and (|20|) with 

\/3 

< fc < — — a, -co < s < oo (22) 

do determine a smooth soliton. This solitary wave is positive, has one maximum 
equal to 1/ (a + 2Va 2 — fc 2 ) and nonzero asymptotes —1/a at x — > ±co, and moves 
with the speed a (a 2 — fc 2 ), that is to the left, as shown in Figures [7] and When 
fc — s- 0, the height of the soliton tends to —1/a, which is the level of the asymptotes, 
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Figure 4. The moving soliton ([19])-([2l]) with a = 1 and k = 0.95: 
t = (solid) and t = 100 (dashed). 
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Figure 5. The shape of the soliton ([1^])— (EU with a = 1 and 
t = 0: k = 0.0001 (solid), k = 0.97 (dot-dashed), k = 0.9999 
(dashed), and k = 1 - 10~ 8 (dotted). 
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Figure 6. The shape of the soliton (jT9j)— ([21]) at < = 0: a = 1 and 
k = 0.0001a (solid); a = 2.0186 and k = 0.97a (dashed). 
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Figure 7. The shape of the soliton, determined by ((19} . ((20} and 
(p?2"} , with a = -1 and i = 0: fc = 0.8 (solid), fc = 0.7 (dashed), 
and fc = 0.6 (dotted). 
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Figure 8. The moving soliton, determined by (JTSJ), (HOI) and ((25}, 
with a = — 1 and fc = 0.75: t — (solid) and t = 4.5 (dashed). 



and the speed of the soliton tends to a 3 . When fc — ► — V3a/2, the soliton's height 
and speed tend to +oo and a 3 /4, respectively. 

The case of "— " sign in the mKdV equation's solution (|15p is processed in the 
same way as the case of "+" sign. As the result, we obtain the following expressions: 



m(x, t) 



a + V a 2 — fc 2 cosh(fcs) 



2fc 2 



ay/a 2 — k 2 cosh(fcs) ' 



(a 2 - fc 2 ) t + as - 4 tanh 1 ( - - Va 2 - fc 2 ) tanh 



fc.s 



(23) 
(24) 



In order to get real, smooth and global solutions of the Qiao equation (QJ , we must 
supplement these expressions either with 



n i ^ 

U < fc < a, — oo < s < oo 

2 



(25) 
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FIGURE 9. The moving antisoliton (|2"3" ]) -(|55 ]l with a = 1 and k = 
0.75: t = (solid) and t = 4.5 (dashed). 




FIGURE 10. The moving antisoliton, determined by ([23]) . (|24|) and 
(|2"6l) , with a = -1 and k = 0.95: i = (solid) and t = 100 (dashed). 



for positive values of a, or with 

\s\ < \ tanrT 1 < fc < -a (26) 

fc -a + Va 2 - fc 2 

for negative values of a. The obtained solitary waves are negative. The antisoliton 
(f25)) - (|2"5"l) moves to the right and has nonzero asymptotes, as shown in Figure 
It corresponds to the soliton, determined by (|T^|) . (p?D]) and (|2"2"|) , in the sense of 
the antisymmetry between positive and negative solutions, discussed in Section [3] 
The antisoliton, determined by (|23|) . (f24)) and (l26l) . moves to the left and has zero 
asymptotes, as shown in Figure [TUJ It corresponds to the soliton (fl"9"|) - (f2"Tj) in the 
sense of that antisymmetry. Owing to this correspondence, there is no need to 
discuss these antisolitons in more detail. 

To conclude the section, we remind the following important observation. When 
we used the soliton solution (fl~5|) of the mKdV equation to obtain smooth real solu- 
tions of the Qiao equation, we had to impose certain constraints on the quantities 
y, t, a and k in the original solution, otherwise the obtained solutions would be 
singular or complex. 
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5. Conclusion 

In this paper, we found a transformation which relates a new integrable nonlinear 
evolution equation of Qiao with the well-known mKdV equation. By means of that 
transformation, we derived smooth soliton solutions of the Qiao equation from 
the known rational and soliton solutions of the mKdV equation. We obtained three 
types of positive solitary wave solutions, or solitons. Solitons of two types have zero 
asymptotes and move to the right, whereas solitons of the third type have nonzero 
asymptotes and move to the left. Also we obtained three types of negative solitary 
wave solutions, or antisolitons, which have the same shape as the corresponding 
solitons but differ from them in their sign and direction of motion. 

It is remarkable that the rational solution of the mKdV equation is transformed 
into a smooth soliton solution of the Qiao equation. An important detail of this is 
that only a part of the rational solution of the mKdV equation produces the entire 
smooth soliton solution of the Qiao equation. When the soliton solution of the 
mKdV equation is transformed, it is also necessary to impose certain constraints 
on its variables and constants, in order that the obtained solution of the Qiao 
equation be real and smooth. These observations may be important for further 
development of the method to obtain parametric representations for solutions of 
nonlinear wave equations via transformations of the equations. 

Last but not least, we believe that the discovered transformation will be useful 
to obtain and study multisoliton solutions of the Qiao equation with the help of 
the Hirota bilinear representation of the mKdV equation [10] . 
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